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Abstract In this paper, we consider initial-boundary value problem of viscoelastic wave 
equation with a delay term in the interior feedback. Namely, we study the following equation 



together with initial-boundary conditions of Dirichlet type in Q x (0, +00), and prove that for 
arbitrary real numbers [L\ and [ii , the above mentioned problem has a unique global solution 
under suitable assumptions on the kernel g. This improve the results of the previous literature 
such as [6] and [13] by removing the restriction imposed on /ii and fi2- Furthermore, we also 
get an exponential decay results for the energy of the concerned problem in the case fix = 
which solves an open problem proposed by M. Kirane and B. Said-Houari in [13]. 
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1 Introduction 

It is well-known that the free vibration of membrane can be described by wave equation of 
the form 



in f2 x (0, +00), subjected to initial conditions and some boundary conditions. Here, £1 C R 2 
is a bounded domain. For the case of high dimension, there are also many discussions in 
mathematical physics. For the better comprehension of our motivation, we appeal to readers to 
keep in mind that the system (1.1) is conservative. 
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Uu{x, t) — Au(x, t) + / g(t — s)Au(x, s)ds + fiiut(x, t) + fj,2Ut(x, t — r) = 
Jo 






To study the propagation mechanism of wave, two main factors, that is damping and time- 
delay, are taken into consideration. When the damping occurs, the corresponding mathematical 
model becomes 

utt(x,t) - Au(x,t) + fxut = 0, (1.2) 

and the energy of this system is dissipative. We here point out that the results about the models 
(1.1) and (1.2) are quite abundant. And what is regretful is that we cannot list these literatures 
one by one. 

Moreover, time-delay effects often appears in many daily life practical problems. These 
hereditary effects are sometime unavoidable, and it may turn a well-behaved system into a wild 
one. That is to say, they can induce some instabilities(see for instance[l-5] and the references 
therein). The stability issue of systems with delay is, therefore, of theoretical and practical 
importance. 

For the wave equation with boundary or interior delay, Nicaise et al. obtained some profound 
results [6-8]. In [6], the authors examined a system of wave equation with the damping and the 
time-delay inside the domain. Namely, they considered the following system 



0, x eft, t > 0, 
x G dft, t > 0, 
x G ft, 

x£fi, < t < r, 



(1.3) 



u tt (x, t) - Au(x, t) + (JLlUt(x, t) + H2Ut{x, t - t) 
u(x, t) = 0, 

u(x, 0) = uo(x), ut(x, 0) = ui(x), 
u t {x,t - t) = f (x,t- r), 

and proved that the energy of the above problem is exponentially decay (or the trivial solution 
is exponentially stable) provided that < fi2 < fJ>i- Furthermore, it is also showed in the case 
^2 > A*i that there exists a sequence of arbitrary small (or large) delays such that instabilities 
occur. The same results were showed when both the damping and the delay act on the boundary. 
Particularly, if fii = 0, that is, if we have only the delay part in the boundary condition, the 
system becomes unstable (see [2]). From these results, we may figure out that the delay term 
lead to instability of the system. 

In this paper, we are concerned with the initial-boundary value problem(IBVP) of the vis- 
coelastic wave equation 

uu(x, t) — Au(x, t) + Jq g(t — s)Au(x, s)ds 

+fiiu t (x, t) + fj, 2 u t (x, t - t) = 0, 
u(x, t) = 0, 

u(x, 0) = uo(x), u t (x, 0) = ui(x), 
u t {x,t- t) = f (x,t- r), 



x efl, t > 0, 
x G dft, t > 0, 
x G fl, 

x eft, o < t < t, 



(1.5) 



where fl C R n ,n > 1 be a regular and bounded domain with a boundary dfl of class C 2 . 
Moreover, fj,±, /j,2 are two real numbers, r > represents the time delay and uo,u±,fo are given 
initial data belonging to suitable spaces. 
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System (1.5) comes from linear models for propagation of viscoelastic wave in a compressible 
fluid. In mechanics, it is well known that solid and fluid materials exhibit not only elasticity but 
also hereditary properties. The hereditary properties are described by viscoelasticity, where the 
mechanical response of the materials is taken to be influenced by the previous behavior of the 
materials themselves. In other words, the viscoelastic materials exhibit memory effects. This 
leads to a constitutive relationship between the stress and strain involving the convolution of the 
strain with a relaxation function. This term is called a "memory term" . From the mathematical 
point of view, these hereditary properties are modeled by integro-differential operators. 

To motivate our present work, let us recall some previous results regarding the viscoelastic 
wave equation. In the absence of time-delay, the viscoelastic wave equation of the form 



in Q x (0,oo), subjected to initial conditions and boundary conditions of Dirichlet type has 
been considered by Cavalcanti et al.[9]. They showed an exponential decay result under some 
restrictions on a(x) and g(t). To be specific, they assumed that a : ft — > R is a nonnegative 
and bounded function, and the kernel g in the memory term decays exponentially. This result 
has been improved later by Berrimi and Messaoudi[10] under weaker conditions on both a and 
g. And then, more general problems than the one considered in [9] are studied in [11,12]. Some 
profound results about the relationship between the stability of solutions and the relaxation 
kernel g are obtained. It is worthwhile to note that, if a(x) = in (1.6), some results about 
the energy decay are obtained by many researchers. See, for instance Tatar [18], Mustafa and 
Messaoudi[19]. 

Then a nature problem is that what would happen when a delay term occurs in (1.6). In 
this case, M. Kirane and B. Said-Houari studied problem (1.5) in [13] with coefficients fi\ and 
[i2 positive. They showed that problem (1.5) has a unique weak solution (see [13], Theorem 
3.1), and its energy is exponentially decay provided that ^2 < Mi- It is worth pointing out 
that the assumptions fj,\ > and ^2 < fJ-i play a decisive role in proof of the above mentioned 
results. Because the authors must use the damping term ix\Ut{x,t) to control the delay term 
(j>2Ut(x,t — t) in the priori estimate of the solution and the decay estimate of the energy. This 
makes the authors arise an open problem that whether the decay property of the energy they 
have obtained are preserved in the case Hi = 0? By the way, the results of [13] has been 
generalized recently by W.J. Liu to a system with time dependent delay in [14]. The method 
used by W.J. Liu is very similar to those used in [13]. 

The aims of the present paper are twofold. At first, we aim to prove an existence result of 
problem (1.5) without restrictions of ^1,^2 > and ^2 < Unlike in [13], our new idea is to 
control the delay term by the derivative of the energy instead of by the damping term in the 
priori estimate of solutions. In the second, we will give a positive answer to the open problem 
proposed by M. Kirane and B. Said-Houari in [13]. That is we will prove a energy decay result 
for problem (1.5) in the case fi\ = 0. The main difficulty in handling this problem is that we have 




(1.6) 
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no damping term to control the delay term in the estimate of the energy decay. To overcome 
this difficulty, our basic idea is to control the delay term by making use of the viscoelasticity 
term. And in order to achieve this goal, a restriction of the size between the parameter fi2 and 
the kernel g and a new Lyapunov functional, different from that one used in [13], are needed. 

The paper is organized as follows. In Section 2 we give some preliminaries. In Section 3, we 
prove the existence and the uniqueness of the solution of problem (1,5). In Section 4, we prove 
the exponential decay of the energy of this problem in the case Hi = 0. 

2 Preliminaries 

In this section, we shall prepare some materials needed in the proof of our result and state. 
We use the standard Lebesgue space L 2 (f2) and Sobolev space -ffd(fi) with their usual scalar 
products and norms. Throughout this paper, C and Cj are used to denote generic positive 
constants. 

We first state the general assumptions on the kernel function g as follows: 
(Al) g : R + — >■ R + is a C 1 function satisfying 

g(0) > 0, 1 - / g(s)ds = I > 0. 
Jo 

(A2) There exists a positive constant £ such that 

g'(t)<-Cg(t), Vt>0. (2.1) 
Nowdet us introduce the following notation: 

(</>oi/>)(t) := f l (p{t-s) f \^j{t) -4>(s)\ 2 dxds. 
Jo Jn 

For any function cf) € C 1 (R) and any ijj G H l (0, T), straightforward computations lead to(see 
[13,15]) 

/ Vt(*) / 4>{t - s)tp(s)dsdx 
Jn Jo 

= l^'o m) - \mmt - ~ {(</> ° m) - (jf ml} ■ w 

The following lemma will be used later in order to define the new modified energy functional 
of problem (1.5). 

Lemma 2.1 (see [13, Lemma 2.2] or [16]). For u G Hq(Q), we have 

J n (J*9(t-*)Wt) -<*))*») dx<(l-l)Cl{goVu)(t), (2.3) 
where C* is the Poincare constant and I is given in (Al). 
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3 Global existence of solutions 



In this section, we will give a sufficient condition that guarantees that the problem (1.5) is 
well-posed. Our main result is the following theorem. 

Theorem 3.1(Global existence). Assume that (Al) and (A2) hold. Then given no £ Hq(Q,), 
u\ <G L 2 (£l), /o € L 2 (£l x (— r, 0)) and T > 0, there exists a unique weak solution u of problem 
(1.5) on (0,T) such that 

u € C ([0,T],H^n)) flC 1 ([0,T],L 2 (O)) , 
Ui € L 2 (o,T;fr^(fi)) nL 2 ((0,T) x ft). 



Proof. Similar to [13], we divide the proof of Theorem 3.1 in two steps: the construction 
of approximations and then thanks to certain prior estimates, we pass to the limit. 
Step 1: Faedo-Galerkin approximation 

Vn > 1, let W n = span{wi, ...,w n } be & Hilbertian basis of the space Hq(Q). We choose two 
sequences (u n o) and (u n i) in W n such that u n o — > uo strongly in Hq(£1) and u n \ — > u\ strongly 
in L 2 {n). 

We define now the approximation 

n 

Un{x,t) = J29jn(t)Wj, (3.1) 
J'=l 

where u n (t) are solutions to the following problem 

In u nU (t)wjdx + J n X7u n (t)wjdx - J* * g(t - s) J n Vu n (s)Vwjdxds 
+ui J Q u nt (x, t)wjdx + u 2 Jq u nt (x, t - r)wjdx = 0, 

U n (x,0) = U n0 (x), U nt (x,0) = U nl (x), 
U nt (x,t) = f n0 (t), t G [-T, 0]. 



(3.2) 



According to the standard theory of ordinary differential equations, the finite dimensional prob- 
lem (3.2) has solution (gj n (t)),j = l,—,n defined on [0,i n ). The prior estimates that follow 
imply that in fact t n = T. 

Step 2: Prior estimates 

Multiplying the first equation in (3.2) by g'j n (t), straightforward computations lead to 



~dt2 



\Unt\\l + ( 1 



J g{s)ds^j \\Vu n \\l + {goVu n ){t) 

= \ ({9 °Vu n )(t) - g(t)\\Vu n (t)\\ 2 2 ) - m J^u 2 nt {t)dx- LL 2 J^u nt {t-T)u nt {t)dx. (3.3) 
By the assumptions (Al) and (A2) about g(t), we have 

(g' o Vu n )(t) - g(t)\\Vu n (t)\\ 2 < 0. (3.4) 
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So, thanks to the Cauchy-Schwartz inequality, we deduce 

J t l (\Wntf 2 + (l ~ J g(s)ds^J \\Vu n \\ 2 + (go Vu n )(t) 



<^J n <t(t-r)dx+(^-^ J a v^{t)dx. (3.5) 



Let 



£(t) = \ (iKtlll + (l - f g(s)ds^j \\Vun\H + {g o Vu n )(t^j . (3.6) 
By integrating on [0, t], we have 

Sit) < jf* jf u l s{s _ r)a!a;( i s + _ ^ jf* u 2 ns (s)dxds + 5(0), (3.7) 

where £(0) is a positive constant. 

Now, using the past history values about u n t(t), t € [— r, 0], the first term in the right-hand 
side of (3.7) can be rewritten as follows 

/ / u 2 ns {s - r)dxds 
Jo Jn 

= In I u np(p) d P dx 

r rO r rt—T 

= y_ u 2 np {p)dpdx + j u 2 np {p)dpdx 

r f0 p rt — T 

= JnJ ^°^ dpdx + J n J u n P (p) d P dx 

<C + [ f u 2 np (p)dpdx, (3.8) 
Jn Jo 



where Co is a positive constant. 
From (3.7) and (3.8), we obtain 



£(t)<(\»2\-»i)J n J o ul p (p)dpdx+ l -^A + £(0)<(\p 2 \-p 1 ) jf S(s)ds + C 1 , (3.9) 

where C\ is also a positive constant. 

Now, we will discuss the inequality (3.9) from the following three cases: 
Case 1: if p\ = \p2\ = P > 0, then Vi > 0, we have £{t) < C\. 
Case 2: if \pi\ < Pi, then Vi > 0, we have £{t) < C\. 

Case 3: if \p2\ > pi, by the Gronwall inequality, once T > be given, then Vi <G [0,T], we 
have S(t) < C(T). 

So, Vpi,p2 € once T > be given, there always exists a positive constant C = 
max{Ci, C(T)} such that 

£(t)<C, Vt€[0,T]. (3.10) 
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Therefore, using the assumption (Al), the estimate (3.10) together with (3.6) give us, for all 
n € N, t n = T, we deduce (u n ) nG N is bounded in L°°(0, T; Hq (fi)), and (u nt ) n ^N is bounded 
in L°°(0,T;L 2 (O)). 

Consequently, we may conclude that u n — u weak* in L°°(0, T; i?Q (Jl)), and « n t — 1 
weak* in L°°(0, T; L 2 (r2)). And we easily deduce that (u n ) ne jy and (u n t) n £N are bounded in 
L 2 (0,T;fl3(fi)) and L 2 (0, T; L 2 (^)) respectively. 

Since the embedding // 1 (0, T; i7 1 (fi)) L 2 (0, T; L 2 (Q)) is compact, we can extract a sub- 
sequence (uk)keN of (u n ) ne N such that Uk —> u strongly in L 2 (0, T; L 2 ($7)). Therefore, uu u 
strongly and a.e. on (0, T) x $7. 

Next, let n = k in (3.2) and pass to the limit, the proof can be completed arguing as in [17, 
Theorem 3.1].D 



4 Exponential energy decay in the case fi\ = 

In this section, we will give the exponential energy decay result about the problem (1.5) in 
the case Hi = 0. That is to say, we will study the following system 

utt(x, t) — Au(x, t) + Jq g(t — s)Au(x, s)ds + fiut(x, t — r) = 0, x € f2, t > 0, 

u(x,t) = 0, x e on, t > o, 

u(x, 0) = uo(x), nt(x, 0) = tii(x), x£fi, 
ut(x,t — t) = fo(x,t — r), x € $7, < £ < r, 

where /u is a real number. 

We first define the classical energy by 



3(*) = ^(lhlll + l|Vu|||), t>0. 



2 

Then by the first equation in (4.1) it is easy to see that 

e'(t)= / Vu t \ g(t — s)Vu(s)dsdx — \i \ u t (t)u t (t — r)dx, t>0. 
Jn Jo Jn 

By (2.2) and inspired by [13], we define the new modified energy functional of problem (4.1) as 
follows: 

E{t) = \ (lM2+ (l-J ||Vn|| 2 + ( 5 oVn n )(t)) + | jf * J^^u^x, s)dxds, 

(4.2) 

where a and £ are two positive constants to be determined later. 

Since the function g is positive, continuous and g(0) > 0, then for any t > to > 0, we have 

f g(s)ds > I' g(s)ds = g . (4.3) 

J0 J0 
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Our stability result reads as follows: 

Theorem 4.1 (Exponential energy decay) Let u be the solution of (4.1). Assume that g 
satisfies (Al) and (A2), and \/j,\ < a. Then, for any to > 0, there exists two positive constants, 
K and A, such that the energy of problem (4.1) satisfies 

E(t) < Ke' x ^' t0 \ Vt > t , (4.4) 

where a, which is only dependant on go, is a positive constant defined by (4.31). 
Proof. We define the Lyapunov functional 

L(t):=E(t) + £l y(t) + e 2X (t), (4.5) 

where £i,i = 1,2 are two positive real numbers which will be chosen later, and 

tf(t) := / uu t dx, (4.6) 
Jn 

X {t) :=- f u t f g(t - s)(u{t) - u{s))dsdx. (4.7) 
Jn Jo 

Remark 4-2 we deduce that, for £; L > 0,i = 1,2, the Lyapunov functional L(t) and the 
energy E(t) are equivalent in the sense that there exist two positive constants /3i,/02 depending 
on £\ and e 2 such that 

fo.E(t) < L(t) < p 2 E{t), Vt > 0. (4.8) 

In fact, using Young's inequality, Poincare's inequality and (2.3), we have £2|x(t)| < ^ll^lli + 
e2(i-i)C'i ^ o and £l \xz,(t)\ < %\\ut\\l + ^UVulll. So, we obtain \L(t)\ < CE(t), where 

C > is a constant. And then (4.8) holds. 

Now, we will estimate the derivative of L(t) according to the following steps. 

Step 1: Estimate of the derivative of E(t). 

Differentiating (4.2), we have 

E'{t) = (u t u tt /^( s ) ds ) VuV ^ - ^)|Vu| 2 ) dx 

[ g(t - s) [ Vu t (t) (Vu(t) - Vu(s)) dxds 
lo Jn 

- I g'(t-s) f \Vu(t) -Vu(s)\ 2 dxds 
2 Jo Jn 

[ u 2 (x,t)dx — - [ e~ aT u 2 (x,t — r)dx 
Jn 2 Jn 

~— t I e-^ t - s K 2 t (x,s)dxds. (4.9) 
2 Jt-r Jn 

Using the first equation of (4.1), we obtain 

E'{t) = \{g'o Vu) (t) - ^(t)||Vu||l - n [ u t (t)u t (t - r)dx 



+ 
+ 
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2" " z 2 
And then, by Young's inequalities, we get 



/ u 2 t (x,t-T)dx- ^ / / e- a ^- s K 2 t {x,s)dxds 
Jn 2 Jt-T Jn 



E'(t) < \ (</ o Vn) (£) + (M + |) | W |1 + _ | e — ) jf n?(*, £ 

-^(t)l|V«||l-^ f I e-^u 2 t (x,s)dxds. (V<5>0) 
2 2 Jt_ T Jf2 



5fep Estimate of the derivative of ^(t). 
Differentiating and integrating by parts, we obtain 



= + J u ^Au -J g(t- s)Au(s)ds - un t (£ - r)j dx 

= IKIll + [ Vti / g(t - s) (Vu(s) - Vu(t)) dsdx 
Jn Jo 

+ ( K J^9(s)ds - lj \\Vu\\ 2 -nj u(t)ut(t - r)dx 

< / Vn / g(t- s) (V«(s) - Vn(£)) dsdx 
Jo 

+ ||«t||i - i||Vu||l - n [ u{t)u t {t - r)dx. 
Jn 

By Young's inequality and Lemma 2.1, we get(see [13]) 

f Vn / #(£ - s) (V«(s) - Vn(£)) dsefa: 
Jn Jo 

< 5i||Vu||l + Q%(£-s)|Vn(s)- Vn(£)|cfc^ dx 

<Si\\Vu\\ 2 + ^(goVu)(t). (V<5i>0) 
4oi 

Also, using Young's and Poincare's inequalities, we have 

-/x / u(t)u t (t - r)dx < 8i\\Vu\\% + C(5i) [ u 2 t {t-r)dx. 
Jn Jn 

Combining (4.12)-(4.14) and choosing S\ small enough, the estimate 

*'(t) < -^||V«||1 + Ci jf («?(£) + u 2 (t - r)) + C 2 ( 5 o Vn) (£) 

holds for some positive constants Ci,i = 1,2. 
5£ep 3: Estimate of the derivative of %(£). 
Differentiating (4.7) and integrating by parts, we have 

x '(t) = (l - [ g(s)ds) [ Vn / g(t - s) (Vn(£) - Vn(s)) dsdx 
V Jo / Jn Jo 



/ g(t - s) (Vm(s) - Vu(t)) ds] dx - [ u t [ g'(t - s) (u(t) - u(s)) dsdx 
Uo J Jn Jo 

- f l g(s)ds\\u t \\ 2 2 + fi f u t (t-r) f g{t - s) (u(t) - u(s)) dsdx. (4.16) 
Jo Jn Jo 

Using Young's inequality and Lemma 2.1, we get(see [13,14]) 

- / g(s)ds) I Vu [ g(t - s) (Vu(t) - Vu(s)) dsdx 
Jo J Jn Jo 



<S 2 \\Vu\\ 2 2 + ^-(goVu)(t), (\/5 2 >0) (4.17) 



-! u t t g'(t- s) (u(t) - u(s)) dsdx < 6 2 \\Vu\\ 2 2 - ^ (g' o Vu) (t) (4.18) 
Jn Jo 02 

and 



H I u t (t -t) [ g{t - s) (u(t) - u(s)) dsdx <^(jo Vu) (t) + S 2 f u 2 t (t - r)dx, (4.19) 
Jn Jo 02 Jn 

where Ci,i = 3,4,5 are some positive constants. 
Combining (4.16)-(4.19) and (2.3), we obtain 

X'(t) < (s 2 - J*9(s)ds^ \\u t \\ 2 2 + 5 2 \\Vu\\ 2 + ^(goVu)(t) 

-ir(9'° Vu) (t) + 5 2 [ u 2 (t - r)dx, (4.20) 
c»2 Jn 

where Q, i = 6, 7 are some positive constants. 
Step 4-' Estimate of the derivative of L(t). 

By using (4. 3), (4. 5), (4. 11), (4. 15) and (4.20), a series of computations yields, for t > to, 
L'(t) < (M + | + £l C 1 + e 2 (5 2 - 5o )) Kill + (e 2 5 2 - ^) ||Vu||2 

+ G - £J f) {g ' ° Vu) {t) + { eiCi + e ^ +lJ f-^)L u2t{t - T)dx 

+ (eiC 2 + £ -^j (g o Vu) (t) - y £ f a e- a(t - s) u 2 (s)dxds. (4.21) 
Now, we deduce that the following system of inequalities 

l -§ + | + eiCi + e 2 (S 2 - g ) < 0, 
e 2 8 2 - £ 4 < 0, 

\ - e -i? > 0, (4.22) 
e 1 Ci + e 2( 5 2 + ^- ^ <0, 

0",£,<5 2 ,£l,£2,<7- > 
is solvable only if we add some suitable conditions to \fj,\. 
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In fact, we can find solutions of (4.22) according to the following steps. 

Step 1: For convenience, we first pick 8=1. And then take 82 small enough such that 

5 2 <- 

and 

5iep As long as 82 is fixed, we select £2 small enough such that 

1 e 2 C 7 



2 8 2 

Step 3: Next, we choose £\ satisfies the relation 



> 0. (4.24) 



50 <^< g0 " 2(52 . (4.25) 



8C1 e 2 2Ci 
By (4.23) and (4.25), we get 

£2(50 - <5 2 ) - eiCi > eiCi + £ 2 5 2 > (4.26) 

and 

£2^2 - ^ < 0. (4.27) 

Step 4'- Now, we must ensure that the first inequality and the fourth one in (4.22) hold. 
That is to say, the following system of inequalities 

+ \i < £ 2 (50 - 82) - eid, 
4H + 2 iM> £ i C 'i+ £ 2<5 2 , (4-28) 
H > 0, o - > 0, £ > 
must be solvable. 

Let fei = E2 (go — 82) — E\C\ and £; 2 = E\C\ + £282- Thus, by (4.26), k\ and A; 2 are two 
positive constants depending on go. And then the system (4.28) becomes 

M + £ < 2fc i> 

-H + ^£>2£; 2 , (4.29) 
M > 0, 0- > 0, £ > 0. 

By (4.26), we have fci > fc 2 . Thus, the system of inequalities (4.29) is solvable only if we choose 
<T small enough such that 



fci e CTr 1 ; 



In fact, the domain 



{(£>|)ei* + x[0,+oo)| 2fc 2 e^<^<2fc 1) 0<H< 2(fcl fc ;f T) } (4.30) 

J. — r 6 
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is the solution space of the system (4.29). 

Therefore, to ensure the solvability of the nonlinear programming problem (4.22), we just 
need to add a condition given by 



H < 1 + e .r = : ( 4 - 31 ) 

here, a is only dependant on g$. 

Consequently, there exist two positive constants 71 and 72 such that 

^p-<- llE (t)+ l2 (g°Vu)(t), Vt>t . (4.32) 

Similar to the steps in [13], the remaining part of the proof of inequality (4.4) can be finished. 
For reader's convenience, we here write the details as follows: 
By multiplying (4.32) by (, we arrive at 

< -7i(E(t) + 72C(S o V«)(t), Vt > t . (4.33) 
Recalling (A2), (4.2), (4.11) and the first inequality in (4.29), we get 

(^<-7i(E(t)- 72 (g>oVu)(t) 

< -ji(E(t) - 2 l2 E\t) + 2fci7 2 ||u t ||! 

< -7iC^(*) " 272^' (t) + 4fc l72J E(t), V* > t - (4.34) 
Now, we add a restriction condition on (, that is, we suppose that 

C > (4.35) 
7i 

Then, there exists a positive constant 73 such that 

F'(t) < - l3 (E(t), (4.36) 

where 

F(t) = (L(t) + 2 72 E(t) ~ E(i). (4.37) 
And then, there exists a positive constant 74 such that 

F'W < -l3(E(t) < -^(F(t). (4.38) 

A simple integration of (4.38) over (to,t) leads to 

F(t) < F(0)e- 74C (*-*°\ Vt > t - (4.39) 

A combination of (4.37) and (4.39) leads to (4.4). The proof of Theorem 4.1 is thus completed. □ 
Remark 4-3 We here point out that the restriction condition (4.35) only makes the admis- 
sible space of the function g a little narrow. But it is not difficult to find that the constant ( can 
be extended to the case ( = ((t), only if ((t) is a positive nonincreasing differentiable function 
with lower bound. 
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